Suppression of Standing Spin Waves in Low-Dimensional Ferromagnets 
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We examine the experimental absence of standing spin wave modes in thin magnetic films, by 
means of atomistic spin dynamics simulations. Using Co on Cu(OOl) as a model system, we demon- 
strate that by increasing the number of layers, the "optical" branches predicted from adiabatic 
first-principles calculations are strongly suppressed, in agreement with spin-polarized electron en- 
ergy loss spectroscopy measurements reported in the literature. Our results suggest that a dynamical 
analysis of the Heisenberg model is sufficient in order to capture the strong damping of the standing 
modes. 



Spin-polarized electron energy loss spectroscopy 
(SPEELS) has recently developed to the point of becom- 
ing a powerful method with which to probe spin waves at 
surfaces and in thin film structures. [1, 2] Landmark ex- 
periments have included the measurement of the magnon 
spectrum of ultrathin Co films on Cu(OOl), [ ] and of a 
single monolayer (ML) of Fe on W(llO). [ ] More re- 
cently, an asymmetry in the magnon spectrum of a bi- 
layer of Fe on W(llO) has also been reported [ ] — a 
direct signature of the Dzyaloshinskii-Moriya anisotropic 
exchange interaction present in this system. [ ] The ex- 
perimental accessibility of these properties challenge the- 
oreticians to address a host of issues in nanoscale mag- 
netic structures, from accurately describing the exotic 
ground states that arise from relativistic spin-orbit cou- 
pling effects, [6-8] to a correct treatment of the dynamics 
of their spin motions. 

In this Letter, we address a phenomenon common to 
all thin ferromagnets consisting of more than one mono- 
layer, which concerns the apparent absence of "optical" 
branches in the measured spin wave dispersion spec- 
tra. [3, 5] In a slab geometry, each layer corresponds to 
an additional atom in the unit cell, and an elementary 
spin wave analysis [9] predicts the occurrence of standing 
wave modes. 

We focus on the fee Co/Cu(001) system, since it has 
been the subject of detailed experimental investigations 
for different Co thicknesses. [ , - ] For a system of 
8 ML, Vollmer et al. [ ] have shown that it can be de- 
scribed by a nearest-neighbour Heisenberg model, char- 
acterized by a dispersion huj = SJS{1 — cos{qao)) along 
the (110) direction, where J is the exchange coupling, 
S is the magnitude of the spin per primitive unit cell, q 
is the length of the magnon wave vector and ao = 2.55 
A. The fit of this curve to the measured SPEELS data 
was JS = 15.0 ± 0.1 meV, which compares well with 
the value of JS = 14.7 ± 1.5 meV obtained from neu- 
tron scattering experiments on bulk fee Co at long wave- 
lengths. [14] Measurements performed on systems as thin 
as 2.5 MLs [^ ^ '^] revealed a weak thickness dependence 



of the magnon energies. However, for all the thicknesses 
reported, there was no trace of any spin wave modes other 
than the lowest energy acoustic branch. 

Previous theoretical investigations of Co/Cu(001) thin 
films, carried out using a random phase approximation to 
a description of the spin response of the itinerant electron 
system, [15, 16] have underlined the role of particle hole 
excitations known as Stoner excitations. Broadly speak- 
ing, these are not relevant at low energies, where the den- 
sity of magnon states dominates, but their role becomes 
dominant at higher energies, suppressing the magnons 
through a process known as Landau damping, [iv] By 
including them in their theoretical analysis, Costa et 
al. [15, 16] successfully described the SPEELS measure- 
ments. In particular, they correctly predict the suppres- 
sion of the optical spin wave modes, and argue this occurs 
primarily through the Landau damping mechanism. Our 
results suggests that drastic damping and broadening of 
the magnon modes is in fact a dynamical feature of the 
Heisenberg Hamiltonian itself. 

Our strategy is to combine first-principles calculations 
with Atomistic Spin Dynamics (ASD) simulations. We 
map an itinerant electron system onto an effective Hamil- 
tonian with classical spins, 

n = -^^Jijmi mj ^K^{mi - exf , (1) 

i/j i 

where i and j are atomic indices, is the classical 
atomic moment, Jij is the strength of the exchange in- 
teraction, and K is the strength of the anisotropy field 
along the direction bk- Other relativistic terms arising 
from spin-orbit interaction are negligible in the case of 
Co on Cu. 

According to the Hamiltonian above, each atomic mo- 
ment experiences an effective interaction field = 
~'im^' temporal evolution of the atomic spins at 

finite temperature is governed by Langevin dynamics, 
through coupled stochastic differential equations of the 
Landau-Lifshitz form, [i6-2Ju] 
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■ ex 
= -7miX[Bi+bi(i)]-7-miX{mi x [Bi + bi(t)]}^ 

(2) 

where 7 is the the electron gyromagnetic ratio, and is 
a stochastic magnetic field with a Gaussian distribution, 
the magnitude of which is related to the phenomenolog- 
ical damping parameter a, which eventually brings the 
system to thermal equilibrium. The coupled equations 
of motion (2) can be viewed as describing the precession 
of each moment about an effective interaction field, with 
complexity arising from the fact that, since all moments 
are moving, the effective field is not static. 

We have calculated the exchange parameters by us- 
ing the method due to Liechtenstein, Katsnelson and 
Gubanov, [2i, 2^] in which the parameters are obtained 
from small angle perturbations from the ground state us- 
ing the magnetic force theorem. [ ] The method relies 
on an adiabatic approximation, in which the slow motion 
of the spins is decoupled from the fast motion of the itin- 
erant electrons. We have determined the exchange values 
based on the relaxed geometric structure, which we cal- 
culated using the projector augmented wave method, [24] 
as implemented in the Vienna ab initio simulations pack- 
age. [25] 

Once calculated, the exchange parameters Jij may 
be Fourier transformed to obtain the so-called adiabatic 
magnon spectrum. In the simple case of a single mono- 
layer (corresponding to one atom per cell), the energy of 
a spin wave with respect to a ferromagnetic ground state 
is given by £;(q) = T^j^o ^Oj [exp (iq • Rq^) - 1], where 
Hij is the relative position vector connecting sites i and 
j. From this expression it is straightforward to calcu- 
late the spin wave dispersion co'(q). [2t)] It is important 
to recognise that in systems with more than one atom 
per cell, the spin wave energies are given by the eigenval- 
ues of an n X n matrix, where n is the number of atoms 
per cell. [ ] Consequently, in a thin film ferromagnet 
consisting of n monolayers, one expects to obtain n spin 
wave branches. 

The principal advantage of combining first-principles 
calculations with the ASD approach is that it allows us 
to address the dynamical properties of spin systems at 
finite temperatures. [i9, 27] Two important quantities we 
focus on in particular are the space- and time-displaced 
correlation function, 



C^(r-r',t) = {m^k{t)m^,k{0))-{m^k{t)){m^,k{0)), (3) 



where the angular brackets signify an ensemble average 
and k the cartesian component, and its Fourier Trans- 
form, the dynamical structure factor 
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Figure 1. (Color online) (a) Adiabatic magnon spectra for 
8 ML Co/Cu(001). Note the optical branches, (b) Compar- 
ison between the acoustic branches for systems of different 
thickness with the SPEELS data reported for 8 ML Co on 
Cu(OOl). [ ,] 
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where q and uj are the momentum and energy trans- 
fer, respectively. 6'(q, cj) is the quantity probed in neu- 
tron scattering experiments of bulk systems [28], and can 
analogously be applied to SPEELS measurements. By 
plotting the peak positions of the structure factor along 
particular directions in reciprocal space, the spin wave 
dispersions may be obtained. [19, 27, 29] 

Fig. 1(a) displays the adiabatic magnon spectrum we 
obtained for 8 ML Co/Cu(001). The most notable fea- 
ture is the presence of several branches, one for each Co 
layer present. At the Brillouin zone center we observe 
several low energy modes. A number of these are energet- 
ically close to the acoustic branch, and naively one would 
not expect these to be washed out by Stoner excitations. 
However, only the lowest acoustic branch, or Goldstone 
mode, is detected experimentally. [ ] As Vollmer et al. 
point out, [ ] this indicates the shortcomings of a static 
interpretation of their data in terms of the Heisenberg 
model. 

Putting these reservations aside for the time being. 
Fig. 1(b) compares the acoustic branches we have cal- 
culated for 1, 2, 3 and 8 ML films with the observed 
SPEELS data. By measuring the curvature of the dis- 
persions of the lowest mode as q ^ 0, we obtain the spin 
wave stiffness, D. The comparison with experiment is 
reasonable: the values of D we obtain are of the order 
of 430 meV A^ - a 15 % overestimate of the experimen- 
tally determined value of 360 meV A^, but considerably 
softer than the theoretical value for 1 ML determined 
by Pajda et al. of 532 ± 9 meV A^, [30] also obtained 
by a real-space adiabatic approach. We also find that 
the magnon energy at the X point in the Brillouin zone 
boundary varies non-monotonically with film thickness 
within the range 260—400 meV. For 8 ML, we find a 




Figure 2. (Color online) Dynamical structure factor *S'(q, ou) 
obtained from ASD simulations of 8 ML Co/Cu(001). Top 
left: Magnon spectrum calculated at 1 K, with a = 3 • 10~^. 
Bottom left: Magnon spectrum calculated at 300 K, with a = 
0.05. Right: Profile of 5'(q, cj), in which the peaks associated 
with the optical branches are more easily visible. 



magnon energy of 300 meV, which is comparable to the 
experimental value of 240 meV. [ ] 

Next, we examine the dynamical behaviour by follow- 
ing the structure factor 6'^(q, cj) defined by Eq. (4), un- 
der different conditions. The top left hand plot in Fig. 2 
displays the spin wave spectrum that results from a cal- 
culation performed at 1 K, in which the damping con- 
stant a in Eq. (2) is set to 3 • 10~^, a value that ensures 
a very weak coupling to the temperature bath. Even in 
this limiting case, in which the role of the temperature 
is minimized, there is little trace of the optical branches, 
especially at small wave vectors. This is made clear in the 
right hand plot of Fig. 2, in which 6'(q, cj) is plotted at 
selected values of q. The bottom left-hand plot displays 
the spin wave spectrum we obtain at 300 K, and using 
a physically more plausible damping constant a = 0.05. 
The contrast with the spectrum directly obtained from 
first-principles shown in Fig. 1(a) is stark: for values of 
q below 0.6, there is no trace of the optical branches. 
They do remain visible at shorter wavelengths, in partic- 
ular when approaching the X point. We also note that at 
300 K the magnon energies at the zone boundary are re- 
duced by roughly 25 meV, relative to the data calculated 
at 1 K, reflecting the influence of temperature effects. 

The suppression of the optical branches near the Bril- 
louin zone center cannot be explained on thermal grounds 
alone, and is a reflection of the dynamical properties of 
the Heisenberg model. To shed light on this issue, we 
consider an atomic chain of spins magnetized along the 
z axis. In this case, the magnon excitations occur only 
along the chain direction {x direction). Fig. 3 depicts the 
eigenmodes projected on the xy plane for three separate 
cases: long wavelength acoustic modes, acoustic modes 
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Figure 3. Magnons in an atomic chain, plotted in the xy 
plane, (a) Long wavelength acoustic magnons; (b) acoustic 
magnons at the BZ boundary; and (c) the optical mode of 
long wavelength magnons for a system with a two atom basis. 
In each case the unit cell is outlined by a dashed box. 

at the Brillouin Zone (BZ) boundary, and the optical 
mode at long wavelengths for a system with a two atom 
basis. The modes shown in Fig. 3(b) and Fig. 3(c) dif- 
fer only in terms of the number of atoms in the unit cell 
and, just as for the Co/Cu(001) system, a calculation of 
the dynamical structure factor for these two modes re- 
sults in a suppression of the optical mode. This effect is 
easily understood by considering the x component of 
in Eqn. (4). When g ^ 0, the sum over cartesian coor- 
dinates for the optical mode (Fig. 3(c)) is exactly zero, 
since it is a sum of alternating parallel and antiparallel 
correlations, leading to a suppression of 5'(q, cj). Note 
that in the case of a one atom basis (Fig. 3(b)), the 
same configuration will instead result in a pronounced 
(and experimentally detectable) structure factor at the 
BZ boundary, due to the phase factor e'^^'^^~^ \ 

This situation is similar to that encountered for 
magnon excitations in rare earth systems with an hep 
structure, [ ] for which it may be shown that the sus- 
ceptibility is 

^(q+r,cj) = i^Ac(q.^)(l+cos(/:)) + i^op(q.^)(l-cos (/:)). 

(5) 

For Co/Cu(001), q is a vector inside the primitive BZ, 
r is an ordering wave vector, and cp = t ■ where p 
is the vector coupling the two magnetic sublattices. If 
q is parallel to the (110) direction, the phase (p is then 
hn -\- kn, where h and k are Miller indices. Inside the 
primitive BZ, the phase (p = and, in the limit q ^ 0, 
Eq. (5) only yields the acoustic term, and it is only this 
mode one expects to detect in experiment. Conversely, 
outside the first BZ (where h = 1 and k = 0) the situation 
is reversed and the optical mode dominates. 

To further illustrate our point. Fig. 4 displays the simu- 
lated spectra of Co/Cu(001) consisting of 2, 3 and 8 MLs. 
In the original experiments the Co films were magnetized 
along the [110] direction, [^] and in our simulations the 
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Figure 4. Dynamical structure factor *S'(q, cj) measured over 
the first and second Brillouin zones, for Co/Cu(001) systems 
consisting of 2, 3 and 8 monolayers. 



precession of the moments is around this axis. The an- 
tisymmetry with respect to the X point is clear, and we 
conclude that the suppression of the optical modes in the 
first BZ is a general phenomenon, in agreement with ex- 
periment. We also note that in a SPEELS experiment, 
the scattered electrons do not penetrate beyond the first 
few monolayers at the surface of the sample. [ ] Simula- 
tions of the structure factor for only the top layer in the 8 
ML Co/Cu(001) structure, reflecting a possible dynam- 
ical measurement of only the top layers, do not differ 
substantially from our results calculated for the entire 
system. 

To summarize, we have studied the experimentally ob- 
served suppression of standing spin wave modes in thin 
magnetic structures, through a combination of first prin- 
ciples calculations and atomistic spin dynamics simula- 
tions. We find that a dynamical treatment of the Heisen- 
berg model is sufficient to capture a drastic dampening 
of these modes near the Brillouin zone center. Addi- 
tional dissipation mechanisms, most notably in the form 
of Stoner excitations at higher values of q, [ 5, 16] ex- 
acerbate the suppression to the extent that in practice 
all trace of these modes is lost in experiment. Given the 
rich spin dynamics of bulk systems characterized by non- 
collinear ground states, [ , ] we suggest their thin film 
equivalents as a promising area for future study. 

We gratefully acknowledge the European Research 
Council (ERG), the Swedish Research Gouncil (VR), and 
the Knut and Alice Wallenberg Foundation for finan- 
cial support. AB and LB also acknowledge eSSENGE 
and SERG, respectively. Gomputer simulations were per- 
formed at the Swedish National Supercomputer Gentre. 
We thank B. Hjorvarsson, J. Kirschner, L. Nordstrom, 
E. Papaioannou and Kh. Zakeri for valuable discussions. 



* andrea.taroni@physics.uu.se 

[1] M. Plihal, D. L. Mills, and J. Kirschner, Phys. Rev. Lett. 
82, 2579 (1999). 

[2] M. Etzkorn, P. S. Anil-Kumar, and J. Kirschner, "Hand- 
book of Magnetism and Advanced Magnetic Materials," 



(Wiley, New York, 2007) Chap. 7. 
[3] R. Vollmer, M. Etzkorn, P. S. Anil-Kumar, H. Ibach, and 

J. Kirschner, Phys. Rev. Lett. 91, 147201 (2003). 
[4] J. Prokop, W. X. Tang, Y. Zhang, L Tudosa, T. R. F. 

Peixoto, Kh. Zakeri, and J. Kirschner, Phys. Rev. Lett. 

102, 177206 (2009). 
[5] Kh. Zakeri, Y. Zhang, J. Prokop, T.-H. Chuang, N. Sakr, 

W. X. Tang, and J. Kirschner, Phys. Rev. Lett. 104, 

137203 (2010). 

[6] L. Udvardi and L. Szunyogh, Phys. Rev. Lett. 102, 
207204 (2009). 

[7] M. Bode, M. Heide, K. von Bergmann, P. Ferriani, 
S. Heinze, G. Bihlmayer, A. Kubetzka, O. Pietzch, 
S. Blugel, and R. Wiesendanger, Nature (London) 447, 
190 (2007). 

[8] M. Heide, G. Bihlmayer, and S. Blugel, Phys. Rev. B 

78, 140403 (2008). 
[9] J. Van Kranendonk and J. H. Van Vleck, Rev. Mod. 

Phys. 30, 23 (1958). 
[10] U. Gradmann, "Handbook of Magnetic Materials," (El- 
sevier, Amsterdam, 1993) Chap. 1. 
[11] C. A. F. Vaz, J. A. C. Bland, and G. Lauhoff, Rep. Prog. 

Phys. 71, 056501 (2008). 
[12] M. Etzkorn, P. S. Anil-Kumar, R. Vollmer, H. Ibach, and 

J. Kirschner, Surface Science 566-568, 241 (2004). 
[13] R. Vollmer, M. Etzkorn, P. S. Anil-Kumar, H. Ibach, and 

J. Kirschner, J. Appl. Phys. 95, 7435 (2004). 
[14] R. N. Sinclair and B. N. Brockhouse, Phys. Rev. 120, 

1638 (1960). 

[15] A. T. Gosta, R. B. Muniz, and D. L. Mills, Phys. Rev. 

B 69, 064413 (2004). 
[16] A. T. Gosta, R. B. Muniz, and D. L. Mills, Phys. Rev. 

B 70, 054406 (2004). 
[17] P. Mohn, Magnetism in the Solid State (Springer, Berlin, 

2003). 

[18] V. P. Antropov, M. 1. Katsnelson, B. N. Harmon, M. van 
Schilfgaarde, and D. Kusnezov, Phys. Rev. B 54, 1019 
(1996). 

[19] B. Skubic, J. Hellsvik, L. Nordstrom, and O. Eriksson, 
J. Phys.: Gondens. Matter 20, 315203 (2008), http:// 
www .fysik.uu.se/ cmt/ asd/. 

[20] J. H. Mentink, M. V. Tretyakov, A. Fasolino, M. I. Kat- 
snelson, and Th. Rasing, J. Phys.: Gondens. Matter 22, 
176001 (2010). 

[21] A. 1. Liechtenstein, M. I. Katsnelson, and V. A. 
Gubanov, J. Phys. F: Metal Physics 14, L125 (1984). 

[22] A. Liechtenstein, M. Katsnelson, V. Antropov, and 
V. Gubanov, J. Magn. Magn. Mater. 54-57, 965 (1986). 

[23] A. Macintosh and O. K. Andersen, "Electrons at the 
Fermi Surface," (Gambridge University Press, Gam- 
bridge, 1980) Ghap. 5. 

[24] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

[25] G. Kresse and J. Furthmiiller, Gomput. Mater. Sci. 6, 15 
(1996). 

[26] J. Kiibler, Theory of Itinerant Electron Magnetism (Ox- 
ford University Press, Oxford, 2009). 

[27] K. Ghen and D. P. Landau, Phys. Rev. B 49, 3266 (1994). 

[28] S. W. Lovesey, Theory of Neutron Scattering from Con- 
densed Matter (Oxford University Press, Oxford, 1984). 

[29] A. Bergman, A. Taroni, L. Bergqvist, J. Hellsvik, 
B. Hjorvarsson, and O. Eriksson, Phys. Rev. B 81, 
144416 (2010). 

[30] M. Pajda, J. Kudrnovsky, I. Turek, V. Drchal, and 
P. Bruno, Phys. Rev. Lett. 85, 5424 (2000). 



5 



[31] J. Jensen and A. R. Mackintosh, Rare Earth Magnetism [33] J. T. Haraldsen, F. Ye, R. S. Fishman, J. A. Fernandez- 

(Oxford University Press, Oxford, 1991). Baca, Y. Yamaguchi, K. Kimura, and T. Kimura, Phys. 

[32] J. T. Haraldsen and R. S. Fishman, Phys. Rev. B 82, Rev. B 82, 020404(R) (2010). 
144441 (2010). 



